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Abstract. We report on experimental transition probabilities for thirteen forbidden [Fe II] lines originating from 
three different metastable Fe ii levels. Radiative lifetimes have been measured of two metastable states by applying 
a laser probing technique on a stored ion beam. Branching ratios for the radiative decay channels, i.e. Ml and E2 
transitions, are derived from observed intensity ratios of forbidden lines in astrophysical spectra cind compared 
with theoretical data. The lifetimes and branching ratios axe combined to derive absolute transition probabilities, 
A-values. 

We present the first experimental lifetime values for the two Fe II levels a^Gg /2 and b^Hn /2 and A-values for 13 
forbidden transitions from a^S 5 / 2 , a^Gg /2 and b"^D 7/2 in the optical region. A discrepancy between the measured 
and calculated values of the lifetime for the b^Hn /2 level is discussed in terms of level mixing. 

We have used the code CIVS to calculate transition probabilities of the a^D-a®S transitions. 

We have also studied observational branching ratios for lines from 5 other metastable Fe II levels and compared 
them to calculated values. A consistency in the deviation between calibrated observational intensity ratios and 
theoretical branching ratios for lines in a wider wavelength region supports the use of [Fe II] lines for determination 
of reddening. 
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1. Introduction 

By definition metastable states cannot decay through elec- 
tric dipole (El) radiation and have therefore long radia- 
tive lifetimes (« 1 s). They are, in general, collisionally 
deexcited in laboratory plasmas. However, in dilute astro- 
physical plasmaa with low probabilities for collisions the 
metastable states can decay radiatively in magnetic dipole 
(Ml) or electric quadrupole transitions (E2), emitting so 
called forbidden lines. Such lines can dominate in emission 
line spectra of low density regions and can be the ma- 
jor cooling agent of dilute plasmas. The use of forbidden 
lines in the diagnostics of a plasma requires the knowl- 
edge of the transition probabilities. These are obtained 
from theoretical calculations since experimental data re- 
quire measurements of both the radiative lifetime and the 
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branching fractions. Neither of these measurements is eas- 
ily performed in the laboratory. 

The first identification of forbidden Fe II lines, [Fe II], 
was done in the pioneering work by Merrill (1928) on the 
emission line spectrum of t) Carinae. The identification 
was stimulated by Bowen’s explanation of the nebular 
lines as Ml transitions between states within the ground 
configuration of 0 III (Bowen 1928). The high spatial res- 
olution of the Hubble Space Telescope (HST) made it pos- 
sible to show that the narrow, strong nebular [Fe II] lines 
in Tj Car emanate from compact gas ejecta (Davidson et al. 
1995), close to the central star. Later HST/STIS spectvd. 
(Gull et al. 1999; Gull & Ishibashi 2001) show the distri- 
bution of [Fe II] in the proximity of the star. It is clear that 
[Fe II] provides possibilities to improve the diagnostics of 
the ejecta of this mysterious star. 

The presence of strong [Fe II] lines in the ejecta of 
Tj Car has initiated the present laboratory study, where 
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Fig. 1. Typical lifetime curves (after corrections) for the levels Fig. 2. HST/STIS spectrum of the blob spectrum of 
and b^Hii/ 2 - Eta Caxinae showing the strong lines a®D-a®S. 




we combine measurements of radiative lifetimes of four 
metastable Fe II levels with theoretical and observational 
data for the branching fractions in order to get absolute 
transition probabilities, A-values. As is the situation with 
A- values (or gf- values) of allowed lines, we will only be able 
to provide experimental data for a few forbidden lines. The 
main bulk of the data for astrophysical analyses must be 
supplied by theoretical calculations. It is therefore impor- 
tant to try to get an estimate of the accuracy of the cal- 
culated data by experiments. The radiative lifetime of the 
upper level of a forbidden transition is the atomic param- 
eter that offers a comparison between theory and experi- 
ment. We have in our experiments chosen such metastable 
states of Fe II which are associated with strong nebular 
[Fe II] lines. One level is known to be affected by level 
mixing (Johansson 1978) and will provide a critical test 
of the theoretical calculations. 

In the present paper we describe briefly the lifetime 
measurements in Section 2. A more detailed discussion is 
given in previous papers (Mannervik et al. 1996; Rostohar 
et al. 2001b). In Section 3 we discuss in detail the method 
of getting branching ratios from astrophysical observa- 
tions by using HST/STIS spectra of rj Carinae and its 
ejecta. New calculations using the CIVS code on the a^D- 
a®S transitions are reported in Section 4. The lifetimes of 
four metastable levels of Fe II and transition probabilities 
for 13 [Fe II] lines are presented in Section 5. In Section 6 
we illustrate the effect of level mixing. 

2. Lifetime measurement of metastable Fe II 
levels 

The atomic structure of Fe II is complex due to seven 
valence electrons in one, two or even three open shells. 
The low complex of even parity configurations (3d 4s) 

has 62 metastable levels in the energy range 0-4.7 eV. All 
of these levels are easily populated in stellar atmospheres 
of moderate temperature, which explains the richness in 


Fe II absorption lines in stellar spectra. The levels are 
also populated in dilute astrophysical plasmas, where they 
decay radiatively by forbidden transitions. The potential 
use of the [Fe II] lines in the diagnostics of cLStrophysical 
plasmas is a strong incentive for measuring the radiative 
lifetime of some metastable states. 

The radiative lifetimes of metastable states are typ- 
ically 6-8 orders of magnitude longer than those of the 
states decaying by El radiation. Consequently, a com- 
pletely different technique has to be used. Firstly, the Fe"^ 
ions have to be stored in a low-density environment during 
the measurement to avoid collisional quenching. Secondly, 
the ions have to be stored under controlled conditions for 
a time comparable to the lifetime. Thirdly, among all the 
Fe“^ ions produced and stored, the ones populating the 
metastable state to be measured have to be unambigously 
selected and recognized. 

A suitable method that meets these criteria has 
been developed at CRYRING, the storage ring facil- 
ity at Manne Siegbahn Laboratory (MSL) in Stockholm 
(Abrahamsson et al. 1993). This method, which is called 
the laser probing technique (LPT) (Lidberg et al. 1999; 
Mannervik 2002), has been applied to Fe II (Rostohar 
et al. 2001b). The Fe+ ions were produced from FeCL in 
an ion source, accelerated up to an energy of 40 keV and 
injected into a storage ring. The long radiative lifetimes 
of the metastable levels require that the ring is evacuated 
to a pressure of 10^^^ torr to avoid collisional quenching 
with the rest gas. The particular state to be measured is 
selected by irradiating the iron beam with monochromatic 
light from a cw laser tuned at a frequency matching an al- 
lowed transition to a higher state. The metastable state 
is thus depleted and excited to a state of opposite parity. 
This short-lived state decays in a fast transition and the 
corresponding fluorescence signal is detected. The laser 
probing pulse is systematically applied at different time 
delays after the injection of the iron beam into the stor- 
age ring. Thus, the number of fluorescent photons is moni- 
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Fig. 3. Branching ratios for the metastable levels &^Gg/ 2 , b^Hn /2 and a®S 5 / 2 . Crosses (x) mark the astrophysical values corrected 
for reddening and the (+) mark the theoretical values (Quinet et al. 1996). 


tored as a function of time and a lifetime curve is obtained. 
Certain correction and normalization procedures have to 
be applied, and all such details are given in previous pa- 
pers (e.g. Rostohar et al. 2001a, b). Examples of lifetime 
curves are given in Figure 1. 

3. Measurements of Branching Ratios using 
Astrophysical spectra 

The relevant parameter in astrophysical analyses is the 
transition probability (Aik) which is related to the mea- 
surable quantities, branching fraction (BElk) and radia- 
tive lifetime (ri) by 

Aik = (1) 

where 

BFik = Aik/ ^ Aik and n = Aik^ . (2) 

The upper and lower states are denoted by i and fc, 
respectively. SFlk is obtained by measuring the intensity 
(in photons per second) of the emission line Aik relative to 
the total intensity of all emission lines from level i. 

Lifetimes of metastable states are typically 6-8 orders 
of magnitude longer than of ordinary excited levels, and 
the transition probabilities of the forbidden lines are ac- 
cordingly very small. To avoid collisional quenching of the 
metastable states a very low electron density (typically 10^ 
cm“^) is required in a plasma emitting observable forbid- 
den lines. In general, this also means that a low concen- 
tration of ions is required throughout an extended source. 
Thus, the combination of low density and low transition 
probability makes dilute, extended astrophysical plasmas 
the best environment to produce forbidden lines of signif- 
icant flux, where the low volume density is compensated 
by a large column density. In rare cases, forbidden lines 
have been observed in laboratory sources allowing accu- 
rate wavelength measurements (Eriksson 1965). 

We have considered two ways to get the BF\s and to 
derive the A- values: 

1) Measure the relative intensities of the relevant [Fe II] 


lines in astrophysical spectra 

2) Use theoretical calculations of forbidden transitions in 
Fell. 

We have finally chosen a combination of the two, but be- 
fore describing the procedure, we comment on the two 
different methods. The main argument against using as- 
trophysical spectra is of course that we want to derive 
the atomic parameter values in an independent way and 
use them in the astrophysical analysis. A practical prob- 
lem with the astrophysical spectra is that the intensities 
are affected by interstellar extinction and may not provide 
the “true” branching fractions. Also, a nebular emission 
line may be blended with other lines from the same or 
another ion. This makes the intensity measurements less 
reliable, even if the blending component is known. The 
general problem with method 2) is the difficulty in esti- 
mating reliable uncertainties of calculated values. 

We will adopt the following approach: We construct 
a reddening curve for one particular spectrum of ejecta 
near 77 Car by using theoretical BF:s for a number of [Fe II] 
lines. We correct the observed intensities of the lines under 
investigation and use them to derive branching fractions. 

3.1. Construction of a reddening curve for rj Carinae. 

We have measured the observed intensities (in photons per 
sec.) of 14 [Fe II] lines in the H ST f ST IS spectrum of two 
compact gaseous condensations (called the Weigelt blobs 
B and D (Weigelt & Ebersberger 1986)) of 77 Carinae and 
grouped the lines according to their upper energy level. 
These lines have been observed at three occasions: March 
1998, February 1999 and March 2000. The line ratios from 
the diff“erent observations show consistent results. We used 
the March 1998 data for the analysis since the lines of 
interest are less affected by blending of other lines. The 
observed intensity ratio between pairs of lines from a com- 
mon upper level, i.e. the uncalibrated branching ratio, is 
compared to the theoretical branching ratio. For the lat- 
ter we have used the calculations by Quinet et al. (1996) 
as they give a good agreement with the experimental life- 
times (Rostohar et al. 2001b). By using a standard redden- 
ing curve (Osterbrock 1989) with a wavelength dependent 
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Fig. 4. A- values for transitions from the upper level b"^D 7/2 for 
different values of the reddening constant C. The theoretical 
values (HFR and SST) are from Quinet et al. (1996). 

exponential degrading of the intensity, we can derive the 
amount of reddening for each line pair. If the exponential 
function is a good approximation the different line pairs 
would give the same value of the exponential constant 
C. By averaging over 11 pairs of [Fe II] lines, we derive a 
value of (7 = 1.6±0.4. The variation in C is caused by line 
blending, by uncertainties in the flux measurements and 
in the calculations of transition probabilities. Some of the 
linepairs used may have larger errors according to a note 
in proof on the paper by Quinet et al. (1996). This has 
been taken into consideration when deriving the uncer- 
tainties of the A- values. However, there is some evidence 
that nebular reddening in the Homonculus may not follow 
the standard interstellar reddening curve, 

3.2. Calibrated branching ratios 

We use the reddening curve to derive calibrated intensities 
and branching fractions for the forbidden lines associated 
with the four met aat able Fe II levels, for which we have 
experimental lifetimes. For one of the levels, a^S, the lines 
fall in a very narrow region (4280-4480 A, see Figure 2) 
so we can, to a first approximation, neglect the effect of 
reddening. We compare the BF's derived from the nebular 
intensities with calculated BF's (Quinet et al. 1996) in 
Figure 3, and find an extremely good agreement. 

We have repeated the same procedure for two other 
levels, b^Hxi /2 and a^Gg/ 2 , which also have their main 
branches within 100 and 200 A, respectively. The agree- 
ment is reasonable good, as shown in Figure 3 even if 
one line from a^G 9/2 is blended. The line at A4177 in 
the ?? Car spectrum is a blend of the two Fe II transi- 
tions a^Fg/ 2 -a^G 9/2 and a^D 5 / 2 -z"^D 7 / 2 - This affects the 
branching fractions for all lines from the a^Gg /2 level. To 
estimate the contribution of the allowed transition to the 
observed nebular feature, we have measured the flux from 
other nebular lines from the same upper level, z^D 7 / 2 - 


Orbital Configurations included Symmetry 

ip 3d^4s4p y^P° 

5p 3d®4p, 3d^5p 

6p 3d®np; n = 4, 5, 6 z®D° 

7p 3d®np; n = 4, 5, 6, 7 z®F° 

5s 3d^4s4p, 3d^5s4p y®P^ 

6s Sd^msnp; m, n = 4, 5, 6 y^P° 

4d 3d®4p, 3d^4s4p, 3d^4s4p4d y®P° 

5d 3d^4p, 3d^4s4p, 3d®4p4d, 3d^4p5d 

4f 3d® 4s, 3d®4p4f a®D 


Based on the calculations by Raassen & Uylings (2000) 
and on solar spectrum observations (Thevenin, 1989), we 
estimate the contribution from the a^D 5 / 2 -z^D 7/2 to 
be ^-^20% of the observed emission feature in r] Carinae. 

Based on this good agreement we can now, in princi- 
ple, combine the theoretical BF:s with measured lifetimes 
for these three levels and derive absolute transition prob- 
abilities (see results in Section 4). 

The branches from the fourth level, b'^D 7 / 2 , are spread 
over a large wavelength interval, 3100 - 5600 A, and the ef- 
fect of reddening is not negligible. In Figure 4 we show the 
transition probabilities derived from the stellar intensities, 
which are corrected for reddening at three different values 
of C located around the adopted value of 1.6. Assuming 
that the standard formula for interstellar reddening is a 
good approximation we see that the error introduced by 
the uncertainty in the C value is rather small. A change 
in C of 0.5 introduces a change in the BF of ~10 % for 
the UV lines around 3400 A and ~20 % for the lines at 
5500 A. 

4. Calculations of transitions probabilities 

The present work uses Cl wave functions of the form 

M 

^{J) = Y.ai^i{aiLiSiJ) (3) 

i=l 

generated by the code CIV3 (Hibbert 1975; Hibbert et al. 
1991). For a specific choice of configuration state func- 
tions (CSFs) {$i}, where the angular momentum cou- 
pling scheme is defined by {oj}, the variationally opti- 
mal expansion coefficients {ui} are the eigenvector com- 
ponents of the diagonalized Hamiltonian whose typical el- 
ement is Hij = < >. The corresponding eigen- 

value is an upper bound to the energy of the particular 
state. In particular, if the eigenvalues {Ei} are ordered so 
that El < E 2 < ' ‘ then 

Ei > (4) 

is a consequence of the Hylleraas-Undheim theorem. 

The Hamiltonian used to determine the final wave 
functions consists of the non-relativistic Schrodinger 
Hamiltonian along with the following relativistic opera- 
tors associated with the Breit-Pauli approximation: mass- 



nl 

5s 


6s 


7s 



fjnZ 

Cjnl 

nl 

^jnl 

^jnl 

Cjnl 

nl 

Cjnl 

^jnl 

Cjnl 

0.07552 

1 

19.02548 

4d 

0.63694 

3 

5.15025 

7p 

14.57355 

2 

0.91682 

-0.27401 

2 

9.35791 


-0.92313 

4 

1.90298 


-35.18975 

3 

1.57708 

0.57776 

3 

5.07754 






18.84919 

4 

1.93717 

-1.03497 

4 

2.43824 

5d 

0.44144 

3 

4.73722 


166.94166 

5 

1.93650 

0.97510 

5 

1.29982 


-0.77547 

4 

2.61070 


-366.00070 

6 

2.15513 





1.04299 

5 

1.28493 


200.95182 

7 

2.40962 

0.06441 

1 

18.38583 









-0.16284 

2 

11.31165 

4p 

5.70277 

2 

1.44819 

4f 

1.00000 

4 

2.44380 

4.62987 

3 

2.03315 


-1.88413 

3 

4.04161 





-10.89075 

4 

2.13989 


-3.85906 

4 

2.46303 





7.93749 

5 

2.07416 









-1.98787 

6 

1.49966 

5p 

0.67165 

2 

2.83628 









-0.78204 

3 

4.21964 





0.97432 

1 

0.30972 


5.72815 

4 

0.85910 





-1.45676 

2 

1.62185 


-6.48385 

5 

1.00256 





11.78939 

3 

1.45012 









-28.95904 

4 

1.58546 

6p 

1.19932 

2 

10.90063 





31.20386 

5 

1.64232 


-1.37372 

3 

10.06883 





-15.20386 

6 

1.64979 


7.29654 

4 

1.35418 





1.40099 

7 

1.11660 


-16.45629 

5 

1.35122 









9.70525 

6 

1.40443 






Table 3. Ab initio and fine-tuned results for calculation including all one- and two-electron changes from (3d® 4s, 3d®4s4d) ®D 
and (3d®4s^ 3d®4s5s, 3d^4p^ 3d®4s4d) ®S. 


Transition Ab initio results Fine-tuned results 


Upper level 

Lower level 

AE (cm-') 

ft 

At (s-') 

AE (cm-') 

ft 

At (s-') 


a®D9/2 

22020.67 

0.0846 

0.9777 

23317.40 

0.0894 

1.6247 


a®D7/2 

21643.88 

0.0837 

0.9030 

22932.64 

0.0886 

1.2042 


3'^Ds/2 

21355.16 

0.0830 

0.6367 

22649.56 

0.0879 

0.8533 


a^D3/2 

21151.27 

0.0825 

0.4061 

22455.02 

0.0875 

0.5469 


a^Di/2 

21029.86 

0.0823 

0.1977 

22340.36 

0.0873 

0.2671 


Table 4. Lifetimes (in s) of metastable Fe II levels. 




b"^D7/2 

a®Ss/2 

o 

to 

b"Hn/2 

Garstang (1962) 
Nussbaumer et al. (1981) 


0.618 

0.326 

0.235 

0.856 

10.1 

Quinet, et ai. (1996) 

SST 

0.567 

0.262 

0.755 

6.59 


HFR 

0.500 

0.220 

0.694 

5.20 

Rostohar et al. (2001b) 
This work (exp) 

This work (calc) 


0.53(3) 

0.23(3) 

0.222 

0.65(2) 

3.8(3) 


SST=Superstructure; HFR=Relativistic Hartree-Fock 


correction, Darwin, spin-spin, spin-other-orbit and spin- 
orbit terms. However, the spin-spin-contact and the orbit- 
orbit terms, which do not contribute directly to fine- 
structure separations, are not included in the calculations. 

The CSFs of all states are constructed from a common 
set of one-electron orbitals of the form 


lPnl(r)Vri^,<P)XmA^) 


(5) 


where the radial functions in Eqn (5) are expressed in 
analytic form as a linear combination of normalized Slater- 
type orbitals (STOs): 


k 

Pnlif) — (^jnlXjnl (^) 
j=l 


( 6 ) 


and where the STOs take the form 


Xjni(r) 


(21, -nO! 




1/2 


exp(~Cjn/^) 


(7) 


Equations Eqn (4) constitute a set of variational princi- 
ples allowing the optimisation of the radial function pa- 


Hamiltonian matrix, i'he integers {Ijni} are kept tixed but 
the exponents {Cjn/} and the coefficients {cjni} in Eqn (7) 
may be treated as variational parameters to be optimised 
subject to the orthonormality conditions: 


f 

Jo 


Pni(r)Pn'i{r)dr = Sn 


I <n^ < 


( 8 ) 


The optimisation of the radial functions on, for exam- 
ple, the energy of just one state would lead to a serious 
imbalance in the accuracy to which the full set of states 
under consideration would be represented. Hence different 
orbitals are optimised on different eigenvalues, so that all 
the states are of a comparable level of accuracy. In the 
present work, the Is, 2s, 2p, 3s, 3p, 3d, 4s functions were 
taken from the ground state Hartree-Fock functions given 
by Clementi h Roetti (1974): hence they are all optimised 
on the a®D state. In the present calculations we used the 
radial functions previously optimised in our calculations 
of the a^D - transitions (Donnelly & Hibbert 2001). 
The additional orbitals consist of 5s, 6s, 4p, 5p, 6p, 4d, 
5d, 4f. The method of optimisation (involving only the 
non-relativistic Schrodinger Hamiltonian) is displayed in 
Table 1 and the radial function parameters are given in 
Table 2. 

We give in Table 3 both our ah initio results of oscilla- 
tor strengths and the corresponding values obtained after 
the ‘fine-tuning’ of the calculated energy levels. This pro- 
cess (Brage & Hibbert 1989) makes small adjustments to 
the diagonal elements of the Hamiltonian matrix so as to 
bring the calculated eigenvalue differences into agreement 
with the corresponding experimental energy differences. 
This process has proved effective in improving the accu- 
racy of calculated oscillator strengths (Hibbert 1996). The 
adopted A-values are inserted in Table 5 and compared 
with the new experimental and previously calculated val- 
ues. 


5. Results and Discussion 

Radiative lifetimes of metastable states in Fe II were first 
measured for two levels, a^Ss/2 and b'^D7/2) and they were 
reported in Rostohar et al. (2001b). In the present paper 
we report on the lifetime measurement of two more levels, 
a^G9/2 and b^Hn/2- All values are presented in Table 4, 
where we also compare with theoretical calculations by 
Garstang (1962), Nussbaumer et al. (1981), Quinet et al. 
(1996) and our new calculations. The theoretical lifetimes 
are obtained by taking the inverse sum of all A-values (see 
Eqn (1)), which is the parameter calculated by the codes. 
Nussbaumer et al. (1981) made calculations only for the 
system of sextet levels, and the good agreement between 
their value for the a®S term and our experimental value 
implies that LS coupling is a good approximation for the 
low-level system in Fe II. Our new calculations of the a®S 
lifetime are also in very good agreement with the measured 
lifetime. 


the relativistic Hcy:tree-Fock (Hk'Kj method in the (Jowan 
code (Cowan 1981) and the Superstructure (SST) code 
(Eissner et al. 1974; Nussbaumer &: Storey 1978). The SST 
code gives systematically larger lifetimes values than HFR, 
but there is no consistency as to which code gives the bet- 
ter value, judged from the experimental values. The values 
calculated by Garstang (1962) are consistently larger than 
the experimental values. 

The branching fractions derived from the emission line 
spectrum of gas ejecta in 77 Carinae have been corrected 
for reddening (see Sec. 3.1). For lines within a limited 
spectral region this correction is very small. For the other 
lines separated by about 2500 A the error introduced is 
the main contribution to the total uncertainty. 

From some of the metastable levels there are a few 
transitions with significant branching fractions that are 
not observed in the nebular spectra. Most fall in wave- 
length regions not covered by the observations. A few are 
weak transitions from the levels of interest that are too 
faint to be observed in the stellar spectrum. For all these 
lines we have used the calculations by Quinet et al. (1996) 
to estimate the missing branching fractions, the sum of 
which results in the residual. The residual is 3.5% for 
a^Gg/2 and 0.1% for a^Ss/2. 

For b^Hii/2 we observe only the two transitions down 
to a^F. There are two more lines with significant branch- 
ing ratios, resulting in a residual of 26%. It is interesting to 
note that the major decay channels (74%) for b^Hn/2 are 
the transitions to a^F, the same channels as for a^G. This 
is a sign of the level mixing between b^H and a^G, which 
is further discussed and illustrated in Section 6. The lines 
responsible for the residual BF appear at 9682 A, which 
is heavily blended, and at 2.03 ju, which is outside the ob- 
served region. The large residual makes it too uncertain to 
derive experimental A-values for the two transitions ob- 
served from b^Hii/2- However, the calculated branching 
ratio for these two lines is in good agreement with the 
observed intensity (see Figure 3). In Table 6 we give ad- 
justed theoretical A-values, i.e. theoretical branching frac- 
tions combined with the measured lifetime. 

As the forbidden lines from b^D7/2 appear in a wider 
wavelength range, 3100-5600 A, the branching fractions 
for the transitions depend on the reliability of the red- 
dening curve discussed in Sec. 3.1. However, instead of 
adopting the calculated A-values for these lines we have 
used the theoretical A-values for 14 lines to derive the fi- 
nal branching fractions from the observed intensity ratios 
in the rj Carinae spectrum corrected for reddening, as de- 
scribed in Sec. 3. Consequently, the uncertainties of these 
BF:s are larger than for the lines from the three other 
levels. The calculated residual is 11% for b^D7/2- 

For the A-values with astrophysical branching ratios 
and experimental lifetimes we have estimated the uncer- 
tainties, which are included in Table 5. The uncertainties 
of the lifetime and the branching fractions are treated as 
independent. In the uncertainty of the BF effects from in- 


Transition 

Upper level Lower level A (A) 

A : a 

Experimental 

Theoretical 




This work 

This work Q“(SST) 

Q‘’(HFR) 

G" 

N“ 

a^G9/2 

a^Fg/2 

4177 

0.29(0.05) (bl) 


0.184 

0.194 

0.14 




4276 

0.83(0.07) 


0.750 

0.819 

0.65 



a^Fs/2 

4352 

0.36(0.04) 


0.345 

0.380 

0.31 


a®Ss/2 

a^Dg/2 

4287 

1.53(0.22) 

1.625 

1.37 

1.65 

1.12 

1.53 


a^Dr/2 

4359 

1.19(0.21) 

1.204 

1.02 

1.22 

0.82 

1.141 


a^Ds/2 

4413 

0.84(0.13) 

0.853 

0.725 

0.858 

0.58 

0.812 


a®Ds/2 

4452 

0.53(0.08) 

0.547 

0.465 

0.548 

0.37 

0.522 


a*Di/2 

4474 

0.26(0.04) 

0.267 

0.227 

0.267 

0.18 

0.255 

b^Dr/2 

a®Dg^2 ' 

3175 

0.23(0.03) 


0.190 

0.208 

0.22 



a^Fg/2 

3376 

0.96(0.10) 


0.846 

0.981 

0.73 



a^Fr/2 

3440 

0.23(0.03) 


0.289 

0.326 

0.24 



a^P5/2 

5551 

0.18(0.04) 


0.149 

0.176 

0.13 



a^Ps/2 

5613 

0.10(0.03) 


0.0808 

0.0934 

0.073 



^’^Quinet et al. (1996); ^Garstang (1962); ^Nussbaumer et al. (1981) 


Table 6. A- values (s"^) for transitions from b^Hn/ 2 . The 
branching fractions from Quinet et al. (SST, 1996) are com- 
bined with the experimental lifetime from this work. 


Transition 
Upper level 

Lower level 

a(A) 

A- value (s ^ ) 

b^Hu/2 

a^Fg/2 

4114 

0.14 


a'Fr/2 

4211 

0.058 


a^Gg/2 

9682 

0.036 


strument calibration, intensity measurements, calculated 
residual and reddening are included. 

6 . Mixing of the a"^G9/2 and the b^Hn/2 levels. 

The largest deviation between experimental and theoret- 
ical lifetimes occurs for the b^Hn /2 level, where the ex- 
perimental value is more than 25% smaller than the calcu- 
lated value. Faster decays or more decay channels are thus 
available for b^Hn /2 than predicted by theory. The result 
is a verification of the level mixing between b^Hn /2 and 
a^Gii/ 2 , which has been observed earlier in the study of 
permitted lines of Fe II (Johansson 1978). The lifetime of 
the a^Gii /2 level is predicted to be about 10 times shorter 
than the lifetime of b^Hn/ 2 - Unfortunately, due to fre- 
quency limitations in the probing laser we have not been 
able to measure the effect due to mixing on the lifetime of 

An astrophysical ’’proof’ for the level mixing is the ob- 
servation of the b^Hii/ 2 - z®F9/2 transition at 6269.97 A in 
the solar spectrum, which contradicts all statements about 
good LS coupling in the lower part of the energy level di- 
agram of Fe II. The transition is a result of level mixing of 
both the upper and the lower states, i.e. b^Hn /2 is mixed 
with a^Gii /2 z®F 9/2 is mixed with z'^F 9 / 2 * Thus, the 
observed line is an appendage to the strong a"^Gn/ 2 “Z^F 9/2 
transition at 5316.61 A (Multiplet 49). The line was not 


identified in the atlas of the solar spectrum by Moore et al. 
(1966). Since the transition probablilty is very sensitive to 
a small change in the eigenvector composition of an en- 
ergy level, there is probably not a need for a large mixing 
to account for a change in radiative lifetimes that matches 
the observed deviation from the calculated value. 

7 . Conclusion 

Experimental transition probabilities of 13 lines from 
three of the four metastable Fe II levels, for which we 
have measured the radiative lifetimes, have been derived 
according to Eqn(l). The values are inserted in Table 5 
together with values from previous calculations. For the 
fourth level we give 3 adjusted A- values (Table 6). We 
also give an estimate of the uncertainty of the transition 
probabilities, see Table 5, which are in the range 10-30%. 
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